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!***■ ' Abstract 

o 

We study the ground state of baryonic/axial matter at zero temperature chiral- 
symmetry broken phase under a large magnetic field, in the framework of holo- 
graphic QCD by Sakai-Sugimoto. Our study is motivated by a recent proposal of 
chiral magnetic spiral phase that has been argued to be favored against previously 
studied phase of homogeneous distribution of axial/baryonic currents in terms of 
meson super-currents dictated by triangle anomalies in QCD. Our results provide an 
existence proof of chiral magnetic spiral in strong coupling regime via holography, 
at least for large axial chemical potentials, whereas we don't find the phenomenon 
in the case of purely baryonic chemical potential. 
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1 Introduction and motivation 

The phase diagram of QCD with temperature and chemical potential is an interesting 
area of both theoretical and experimental studies. On the experimental side, we have an 
access to some part of regions through heavy-ion collision and astrophysical observations 
such as neutron stars. However, extracting reliable information on the state of QCD 
matter from these observations is not quite trivial due to theoretical difficulties in dealing 
with strong interactions. For asymptotically high temperature or chemical potential, 
we seem to have reliable perturbative QCD computations to predict deconfined gas of 
quark/gluons and color-flavor-locked superconducting phase [1] respectively. However, 
perturbative QCD picture often seems to break down for interesting regimes such as quark- 
matter inside neutron stars and quark-gluon plasma in heavy-ion experiments [2]. In fact, 
these situations get more complicated than the simple (T, /i) phase diagram, because they 
typically involve extremely large magnetic field ranging roughly 10 12 — 10 18 G (gauss), and 
the QCD matter can feature completely new behaviors under these conditions. Examples 
include magnetic catalysis [3], anomalous axial current in dense matter [3], chiral magnetic 
effect [5], new mechanism of pulsar kicks [6], spontaneous magnetization of neutron stars 
[7], etc to name a few. 

As first-principle QCD computations are hard in most cases of interest, it is worthwhile 
to explore other indirect methods to attack the problems, and the AdS/CFT-inspired 
approach of holographic QCD is one such example. There has been recent outburst of 
research in this direction applied to many aspects of QCD. For our more specific subject 
of studying effects from magnetic field, previous works include Ref . [El E| QUI El El EJ 
□3 CEHJ HE, Ej. One interesting feature in the studies in Refs. pU HI H21 El HI US, E] is 
the interplay between magnetic field and triangle anomalies of global symmetries of QCD. 
See also Refs. [T8l [T9l 120] for hydrodynamic effects with triangle anomalies. In all these 
studies, one is interested in physics of magnetic field in the deconfined phase of QCD 
plasma at high temperature. 

The physics of QCD with magnetic field in a low temperature, chiral-symmetry broken 
phasej has also been of much interest in the field theory side. One question that was 
addressed in Ref. [23] is what would happen to the QCD vacuum when one introduces 
small baryonic chemical potential in the presence of magnetic field. Naively, if the baryon 



*We focus on the Sakai-Sugimoto model [5T] where chiral-symmetry restoration happens simultane- 
ously with deconfinement. We postpone the study of more general cases in Ref. [35] to the future. 



chemical potential is much smaller than the nucleon mass of GeV, the real nucleon won't 
play a dominant role in the phase, and the physics of chemical potential might be more or 
less empty. However, Ref. [23] observed that in the presence of magnetic field, the effect 
from axial anomaly can change the picture. Within the low energy effective theory of 
chiral Lagrangian where pion field U(x) = exp ( ■A-n(x)) is a relevant degree of freedom, 
the conserved gauge- invariant baryon number current is given by 

B» = — —e tiVaP Tr(u- 1 d v UU- 1 d a UU- 1 dBU) 
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L ; e^d v hAl Mr Ti (Q(U- l d p U + dpUU- 1 ))] . (1.1) 
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where the first term is the topological Skyrmion number current representing real nucleons 
as Skyrmion solitons, while the second piece of our interest can be traced back to the Wess- 
Zumino-Witten term due to axial anomaly. It should be emphasized that the above result 
is uniquely determined by requiring both conservation and gauge-invariance of the current 
[2"4"] . It is also worth of mentioning that holographic QCD is an ideal set-up to reproduce 
the above result, because 5-dimensional bulk gauge field unifies the dynamical field of 
pions and external gauge potential coupled to the current in a single framework 
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(U- X QU) ^+ + QV-] A^ M + ip+U^dpU + (excited modes) , (1.2) 



so that the holographic version of baryon current which is automatically conserved and 
gauge- invariant 

I r+oo 

B» = — &z<T°* Tr (F va F fiz ) . (1.3) 

must reproduce fll.ip . which is indeed the case [23]. With the second term in ( II. ip . one can 
have a non-zero baryon charge density with a magnetic field B and a gradient of pion 
U~ 1 dU ~ dn along the same direction, without creating heavy nucleons as Skyrmions. In 
the massless QCD with sufficiently small baryonic chemical potential, it has been argued 
that this phase with constant gradient of pions is indeed the favored ground stateJH As 
the phase of chiral condensate U(x) is rotating along the direction of pion gradient, this 
phase is often called chiral spiral phase. Note that this name also appears in several 
1+1 dimensional models of chiral symmetry breaking such as chiral Gross-Neveu model 
(2- dimensional Nambu-Jona-Lasinio model) or t'Hooft model in the presence of chemical 
potential for the same reason [26J. In fact, it is possible to make a connection : for 



^If quarks have non-zero current masses, the pion gradient would not be constant, rather make domain 
walls [22]. 



sufficiently strong magnetic field, quarks are confined to the lowest Landau level and the 
system effectively behaves like 1+1 dimensional system^ 

Interestingly, there recently appears a new proposal by Basar-Dunne-Kharzeev [28] 
regarding the true ground state of the system, at least for sufficiently strong magnetic 
field. Based on the effective 1+1 dimensional reduction of quarks to the lowest Landau 
level and assuming existence of various bi-fermion condensates due to QCD interactions 
consistent with 1+1 dimensional Lorentz symmetry, they argue that one generically has 
condensates of transverse currents to the direction of magnetic field, in addition to the 
above mentioned chiral spirals of pion which in fact corresponds to an axial current 
dictated by axial anomaly. The fermion pairings responsible for these condensates have 
net non-zero momentum, so the condensate of transverse current looks like spirals and 
inhomogeneous. It is important to mention that they considered both baryonic and axial 
chemical potentials, and anomaly dictates existence of axial/baryonic currents along the 
magnetic field respectively. The latter phenomenon is chiral magnetic effect |5]. As 
the transverse currents resemble chiral magnetic current along the magnetic field due to 
triangle anomaly, these spiral-shaped transverse currents were named as chiral magnetic 
spirals. 

Although their arguments seem quite generic, the basic underlying assumption is the 
validity of quark picture which might not be clear in the confined/chiral symmetry broken 
phase we are considering^! Another point is how to deal with strong interactions which 
become non-perturbative in the regime of interest. As holographic QCD is aimed at dual 
description of strongly coupled regime of QCD, it seems interesting to confirm/disprove 
the chiral magnetic spirals within the set-up, and this is the main purpose of this paper. 
Our work may be considered as a strong-coupling test of the proposal. For our objectives, 
it is sufficient to work in the holographic model of Sakai-Sugimoto [2T] as it features all 
the relevant chiral symmetries of QCD as well as its triangle anomalies. 

In the model of holographic QCD by Sakai-Sugimoto, Ref.fTOt [TT] constructed the 
solution which precisely maps to the previous homogeneous chiral spiral phase. See also 



•"■Even without magnetic field, there have been proposals of effective 1+1 dimensional reduction on 
the Fermi surface of dense quark matter leading to quarkyonic chiral spiral phase |27j . We however focus 
on the case of low chemical potential with magnetic held in the confined phase. 

Hn this respect, a very strong magnetic field beyond QCD scale might help because quarks can be 
restricted to a Landau cell whose scale can be short enough for the validity of asymptotic freedom [TU] . 
We are more interested in the regime where this is not necessarily the case, but still would like to ask the 
question of existence of chiral magnetic spirals. 



Ref . [12| for more study on the solution including iso-spin chemical potential as well. 
Axial anomalies in QCD manifest themselves as 5-dimensional Chern-Simons terms in 
the holographic QCD, and indeed they play a crucial role in the existence of the solution. 
An immediate generalization of their solutions is to include axial chemical potential as 
well as baryonic one, and this is also what we would like to consider in searching for chiral 
magnetic spirals. Although axial symmetry is not strictly conserved due to anomaly with 
QCD gluons, the violation is -^-suppressed in large N c limit that holographic models of 
QCD are based on. A more conservative viewpoint is to take axial chemical potential or 
axial charge density as a temporary, dynamical fluctuation localized in space-time. The 
"phase" with axial chemical potential should then be viewed as dynamical/local phase 
in the situation where these fluctuations are sufficiently slowly varying with respect to 

In this more general homogeneous chiral spiral solution of the Sakai-Sugimoto model 
with baryonic/axial chemical potentials and the magnetic field in the chiral-symmetry 
broken phase, one can easily check the triangle-anomaly dictated induced currents along 
the magnetic field. If the proposal of chiral magnetic spiral is valid at least in some 
region of parameter space of chemical potentials and magnetic field, this homogeneous 
chiral spiral solution in fact won't be a stable ground state against instability towards 
forming chiral magnetic spirals, i.e. transverse currents with net momentum. Our idea is 
therefore to search for any linearized instability of transverse currents in the fluctuations 
from the given homogeneous chiral spiral solution. We will first give a general argument 
for the existence of this instability towards chiral magnetic spirals with sufficiently large 
axial chemical potential, and provide numerical evidence of our phase diagram of chiral 
magnetic spirals in the space of baryonic/axial chemical potentials. This is at least an 
existence proof of the phenomenon in strong-coupling. The Chern-Simons term is essential 
to have such an instability, so the phenomenon is one dynamical effect of triangle anomaly. 
We also point out that the case of purely baryonic chemical potential which was the subject 
of previous studies is safely stable against the instability, so our results indicate a closer 
tie between axial chemical potential and chiral magnetic spirals. This deviates slightly 
from the original argument of Ref. (28] where both baryonic and axial chemical potentials 
are causing the phenomenon, and it seems interesting to understand the difference better. 

It is worth of mentioning that several similar instabilities of transverse currents have 
previously been found in holographic models [29j E3 EH [32] and the basic mechanism 
is universal : strong electric field plus Chern-Simons term can induce instability |30j . 



In terms of mechanism, our example also belongs to the same category. However, our 
situation is much more delicate because of co-existence of magnetic field and electric field 
that are correlated with each other by equations of motion of the model. In fact, this 
is why it is physically more interesting due to astrophysical relevance in magnetars or 
spinning neutron stars and off-center heavy-ion collisions. Another realization of this 
aspect is the non-trivial phase diagram of chiral magnetic spirals we present in this work. 

2 Review on magnetized axial/baryonic matters in 
the Sakai-Sugimoto model 

In this section, we review the holographic solutions of magnetized axial/baryonic matters 
in the Sakai-Sugimoto model in Refs. [I~0| ITT] , which will serve as a starting point for our 
search for chiral magnetic spirals of transverse currents to the applied magnetic field. 
To be more precise, we follow Ref. [10J which constructed the solutions of two derivative 
truncation of the full DBI action. More elaborate solutions of Ref. [TT] with the full DBI 
may be used in a future study. Readers who are familiar to these solutions may skip 
details of this section. 

The Sakai-Sugimoto model of holographic QCD is based on pairs of probe D8-D8 
branes for each quark flavor of QCD. These probe branes, meaning that their back- 
reactions to the background geometry are neglected (which is parallel to quenched ap- 
proximation in lattice QCD), are embedded in some background geometry produced by 
5' 1 -compactified N c D4 branes. These D4 branes or the geometry they produce are re- 
sponsible for SU(N C ) gluonic dynamics of QCD. The dynamics of D8 branes corresponds 
to chiral dynamics of flavor quarks. The description becomes simple in large N c and large 
t'Hooft coupling A = g\ M N c limit where one can neglect various quantum and higher 
derivative corrections in this dual picture. The model realizes the chiral symmetry break- 
ing of massless QCD SU(N F ) L x SU(N F ) R — y SU(N F ) V in an intuitive way by adjoining 
D8-D8 branes each representing SU(N F ) L and SU(N F ) R respectively, so that the residual 
symmetry on D8 branes at IR is simply SU(N F )v- One can also study finite temperature 
deconfined/chiral-symmetry-restored phase of the model. Relevant chiral symmetries in- 
cluding their axial anomalies are consistently packaged into 5-dimensional effective action 
on D8 branes world-volumes. 

Because the main features of the physics we are interested in exist even in the case of a 
single flavor N F = 1, we will study a single pair of D8 and D8 branes for simplicity, whose 



extension to larger number of flavors should be straightforward. This wouldn't necessarily 
mean that the case of multiple flavors would have same results as those results we obtain 
for a single flavor, because non-Abelian iso-spin symmetry might introduce a new feature 
that is absent in our analysis. We also simplify our analysis by working with two derivative 
expansion of the full Dirac-Born-Infeld action as in Ref.pl)]. Although this would imply 
that the detailed analysis and results of the present work have some limitation to not- 
too-big magnetic fields, we expect that general conclusions we draw from the present 
analysis won't change much in the full theory of Dirac-Born-Infeld action. We leave these 
questions as open problems for the future. 

The action of D8-D8 system up to two derivatives is [33jj 

S = -kJ dx 4 dZ [i (l + Z 2 ) ~* F, U F^ + M 2 KK (l + Z 2 ) F, Z F» Z 

+ ^- 2 Jdx*dZe MNp Q R A M F NP F QR , (2.4) 

where k — f ( M ) . We are interested in a zero temperature, chiral-symmetry broken 
phase where two D8 branes join at IR to form a single D8 brane. The second line is the 
Chern-Simons term responsible for triangle anomaly of U(1)l x U(1)r broken down to 
U(l)v- Note that U(1)l x U(1)r also suffers from QCD instanton-induced anomaly with 
two external gluons, but it can be ignored in large N c limit in a first approximation. The 
equations of motion are 

1 + Z 2 Y l d^ + M 2 KK d z \(l + Z 2 ) F Zv \ - J^— e ^F Za F M = , 

M 2 KK (l + Z 2 )d^ z + -^^e^F^F a$ = 0, (2.5) 
and our subsequent analysis will be largely based on these only. 

To realize magnetized axial/baryonic matters with an external magnetic field pointing 
to x 3 -direction with a strength B, one looks for a solution with only Fq Z , F 3Z and F\2 
being turned on with an assumption of homogeneity d^ = at the field strength level. 
This is motivated by the corresponding field theory situation with baryonic chemical 
potential studied by Son-Stephanov [23] where the solution has a constant gradient of 
pion along x 3 ( in our simplified case of Np = 1, it would be U(1)a ?7-meson). Because in 
the Sakai-Sugimoto model, the gradient of pion field enters as 

lim Z 2 F 3Z + lim Z 2 F 3Z ] ~ J^ xial ~ d 3 n + (axial vector mesons) , (2.6) 

Z—too Z— >— oo / 



^Greek indices span the Minkowski (3+l)-dimensions, while capital letters will be used for the full 
5-dimensions. Our metric signature is T) = ( — h + + +) and the numerical e symbol is e 0123Z = +1. 
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one naturally expects that the holographic solution would be homogeneous in terms of 
field strengths. Note that the above statement is completely gauge- invariant, while if one 
chooses to work with A 5 = gauge, then d^n will also appear in the boundary values 
of (^(oo) — A^—oo)). However, one should note that this is true only in the specific 
^.5 = gauge, and is not a correct gauge-invariant statement. This clarification becomes 
important when one considers both axial and baryonic chemical potentials, as the presence 
of axial chemical potential will induce contributions from vector mesons too which can be 
read off from the field strengths only, 



lim Z 2 F 3Z - lim Z 2 F 



rz 



j-baryonic 
J 3 



(vector mesons) 



(2.7) 



y Z— 5>oo " Z— ►— 00 

The induced j 3 ar y° mc j n terms of vector mesons due to axial chemical potential can be 
identified as chiral magnetic effect in the chiral-symmetry broken phase we are considering. 
With this Ansatz, the equations of motion simplify as 



dz 
Oz 



l + Z 2 )F 



:vz 



Nr. 



8n 2 nM 2 



-FqzF\2 



KK 



1 + Z 2 )F 



oz 



N r 



8ti 2 kM 2 kk 



FzzF\2 



while one of the Bianchi identities tells us 



(2.1 



dzF-10 — doF-\ 



'2-T1Z 



dxF, 



2Z 



(2.9) 



so that F 12 is a simple constant along the radial direction Z, 

F12 = B (constant) 
Using this fact, the general solutions of (12. 8 j) are easy to be obtained as 



Fqz 

F3Z 



— (c A cosh [stan _1 (Z)l + C B smh fstaa _1 (Z) 

1 + Z 2 v L J L 



1 



1 + Z 2 



C B cosh Btan' L (Z) +C A sinh Btim~ l (Z) 



(2.10) 



(2.11) 



with B 



N C B 



&-k 2 knP — • ^ ne ^ wo integration constants Ca and Cb are free at the level of equa- 
tions of motion, but ultimately they should be determined by two chemical potentials, 
axial \Xa and baryonic [ib respectively. Note that previous literatures considered only 
baryonic chemical potential //# ( or equivalently Cb ), although the above is a straight- 
forward generalization of theirs. (j2.1ip represents a holographic situation where one has 



a magnetized matter with axial/baryonic chemical potentials without having real baryons 
excited. From a non-zero F 3 z, one also notices an induced axial/baryonic current along 
the magnetic field, and especially axial current is represented by a constant gradient of 
pion field. Physically this phase will be favored against exciting real baryons, when bary- 
onic chemical potential is smaller than the nucleon masses ps ^ Mat ~ 1 GeV to create 
real baryons. Recall that one is at zero temperature in a chiral-symmetry broken phase. 
The reason one can have baryonic/axial charges without having real baryons is in fact 
due to anomaly : external magnetic field plus an induced axial/baryonic current carry 
non-zero baryonic/axial charge density, respectively. As this is from triangle anomaly of 
U(1)l x U(1)r, its holographic realization is played by the 5D Chern-Simons term, as one 
can easily observe in the above solution. 

It is not difficult to relate one integration constant Ca with the axial chemical potential 
Pa- From the definition, one has 

1 1 1 /■ 

Pa = --(p L - Pr) = -^(A (+oo) - A (-oo)) = --J dZF z 



-i fdZdz 
2BJ 

sinh (zjs) 
B 



zo 

Z J 

1 + Z 1 ) F 3Z \ 



1 

IB 


lim (l + Z 2 ) F 3Z - 


-lim (i + Z 2 )f 3Z 

Z— >-oo \ / 



C A , (2.12) 



where one has used the equation of motion (J2.8P between the first and second lines. Hence, 
Ca can be thought of as the axial chemical potential up to a constant factor. Interestingly, 
it is much more subtle to identify correct relation between Cb and the baryonic chemical 
potential ps- Note that the solutions presented above do not give ps straightforwardly, 
contrary to the axial chemical potential ; pb as the value of |(Ao(-|-oo) + A (— oo)) can't 
be determined from the solution (12. lip itself. There have been two proposals how to relate 
Cb with pb- In Ref.jlO], they proposed to minimize the free energy (H — PbNb) with 
respect to Cb given a fixed value of ps, to have a relation between Cb and ps- Here Nb 
is computed from the solution (12. lip by the Chern-Simons induced baryon number 

Nb = ^I d " xdz eiJkF v FkZ = lir / Sxdz F ^ > ( 2 - 13 ) 

while the Hamiltonian H is computed from the Yang- Mills part of the action only. Note 
that in this prescription, one needs only field strengths to compute things. Another 
proposal in Refs. [Til [12] is to extremize the on-shell action including the Chern-Simons 
term, instead of the previous (H — pbNb)- Notice that since Chern-Simons term requires 



gauge potential in addition to field strengths, one has to integrate the solution 02. 11 j) 
to get the gauge potential, and the \ib as the value of ^(A Q (+oo) + A (— oo)) enters 
the expressions in this way. Presumably the two prescriptions should have agreed, but 
curiously they give disagreeing results with each other. The basic reason for this ambiguity 
seems to be non gauge-invariance of the 5D Chern- Simons term. A similar kind of issue 
with 5D Chern- Simons term has also been observed recently in computing zero- frequency 
chiral magnetic conductivity in a deconfined phase [15] . See Refs.JITJ [MJ [35j [36] for 
proposals of the resolution. To be conservative, we therefore leave Cb as it is and present 
our results of analysis in terms of (Cb, Ca) parameter space. One can easily translate our 
results in terms of (hb, Ha) instead of (Cb, Ca) after the issue will be settled down in the 
future. It is however always true that (Cb, Ca) are monotonically increasing functions of 
(Hb, Ha), so one can easily get a qualitative picture from our results. 

3 Search for holographic chiral magnetic spirals 

We come to our main motivation of looking for a signal of chiral magnetic spirals proposed 
by Basar-Dunne-Kharzeev [28] . According to the proposal, the homogeneous solution we 
describe in the previous section may not be the most favored state of the system, at 
least in some region of the parameter space (Cb,Ca)- One instead has condensates of 
transverse axial/baryonic currents which possess a finite momentum along the magnetic 
field, so that profiles of these currents form spiral shapes. Their argument is based on an 
effective reduction of the system to (l+l)-dimensions in the presence of strong magnetic 
field with perturbative dispersion relations of quarks/anti-quarks, although formation of 
various bi-fermion condensates is assumed to happen due to QCD interactions. As QCD 
is strongly coupled in the chiral-symmetry broken regime, it would be valuable to check 
the proposal independently in the framework of holographic QCD, which can serve as an 
alternative proof of the existence of the phenomenon in a strongly-coupled regime. 

To see whether chiral magnetic spirals are more favored than the homogeneous phase 
in the previous section, we study linearized stability of the homogeneous solution against 
small perturbations of transverse currents. In case one finds unstable modes of finite mo- 
mentum along the magnetic field, it is a sufficient proof of the existence of chiral magnetic 
spirals. It will also be very interesting to construct the end point of this instability, but 
we envision that this study may be possible only numerically and we leave this task as 
another future direction. 



The homogeneous solutions (12. lip preserve a residual rotation symmetry in the trans- 
verse (x 1 ,x 2 )-plane, so that linearized fluctuations can be classified according to their 
helicity, and different helicity modes are decoupled from each other. The helicity modes 
are fluctuations from (A Q ,A 3 ), and we are not interested in them in the present study. 
Instead we focus on the helicity ±1 transverse modes from (Ai, A2), which will be respon- 
sible for a formation of chiral magnetic spirals. Starting from the equations of motion 
(12. 5p . it is straightforward to obtain the linearized equations of motion governing fluctua- 
tions. Let us look at particular modes with definite frequency and momentum, e - tU]t + % P x ; 
with real p and possibly complex u) to allow instability if it exists. Defining helicity ±1 
modes as 

A& = 8 A x ± iSA 2 , (3.14) 

they satisfy the linearized equations 



(l + Z") _i (u; 2 - p*) #) + Ml K d z [(l + Z 2 ) d z A&\ ± -^f- ( P F 0Z + coF 3Z ) #> = 



(3.15) 



where F oz and F 3Z are background solutions of ( 12.1 If) . This is our master equation that 
we focus in the following. Note that since we have a time-reversal symmetry t — > — t in 
our background (this contrasts to the case of finite temperature with black-hole horizon 
where in-coming boundary conditions break time- reversal symmetry), a solution with uj 
will always accompany another solution with its conjugate Q. Like a harmonic oscillator 
with unstable tachyonic potential, an instability can be identified simply from having a 
complex- valued u) with non-zero imaginary part, so that one mode out of a; or a) grows 
exponentially in time. One also observes that the last term in the equation coming from 
the Chern-Simons term is chiral depending on the sign of helicity, and this term will in 
fact be crucial to have a possible instability. The role of Chern-Simons term in causing 
instability was already noticed some time ago by several authors, initiated by Domokos- 
Harvey [29]. The basic lesson is that in the presence of sufficiently strong electric field, 
Chern-Simons term induces instability, which was nicely clarified by Nakamura-Ooguri- 
Park recently [30] . See also Ref . j3l] for another realization of the instability in the Sakai- 
Sugimoto model in a different set-up from ours, Ref. [32] for a generalization to higher 
dimensions, and Ref. [37] for 2+1 dimensional system for condensed matter application. 

Perhaps, the main difference of our situation from previous cases is that we have both 
electric field Fq Z as well as magnetic fields F 3Z and F12, whose strengths are correlated 
with each other, and in fact, this correlation is precisely why the current situation is 
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physically interesting, such as magnetars or spinning neutron stars with strong magnetic 
fields. One can easily notice that what matters for instability is a Lorenz scalar F 2 in 5D, 
not electric or magnetic fields alone, so that one can get an interesting phase diagram of 
instabilities in the current situation, which is what we aim at in this work. 

Given (Cb, Ca) which enters (13. 15[) as F oz and F 3Z , one studies an eigen-mode problem 
of u for all possible real momentum p with normalizable boundary conditions on A^ at 
the UV boundaries Z — > ±00, 

A {±) ~ - , Z ->■ ±00 . (3.16) 

The chiral magnetic spiral will be present if u is complex for some finite range of values of 
momentum p. In general, one needs numerical analysis to map out the region of (Cb, Ca) 
having chiral magnetic spirals. 

Before presenting the results of our numerical study, it is possible to estimate a rough 
picture how the phase diagram in (Cb, Ca) would look like. One first notices that if one 
removes the last term in (j3.15p . the equation becomes identical to the Kaluza-Klein reduc- 
tion equation of massive vector mesons in the Sakai-Sugimoto model, whose eigenvalues 
A 2 of the operator 



d z [[I + Z 2 )d z A\ =-A 2 (l + Z 2 ) "A , (3.17) 

are mass squares of massive vector mesons in units of M\ K . From previous studies, 
A 2 > X 2 min = 0.669 where \ 2 min is the mass square of the rho meson [21]. As a rough 
approximation, one therefore replaces this operator in (I3.15P by 

M 2 KK d z [(l + Z 2 ) d z A&] ~ -M 2 KK X 2 eff (l + Z 2 ) ~ h A& , (3.18) 

with some \ 2 ^ > \ 2 min - Next, because of the UV boundary conditions A^ — > as 
Z — > ±00, one can assume that the wavefunctions are localized near Z = at least for 
low-lying states, and as a crude approximation one therefore substitutes F oz and F 3Z in 
(13.15P by their values at Z = 0, 



oz 



C A , F ZZ ~-C B ■ (3.19) 



Then the equation (13. 15[) becomes an algebraic equation for ui and p in this rough ap- 
proximation, 

(co 2 -p 2 )-\l ff M 2 KK ±-^-(C A p-CBU>) = , (3.20) 
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No instability 




Figure 1: A rough estimate of the expected phase diagram of chiral magnetic spiral in 
(Cb, Ca) space. 



which can be brought into the form 



WT 



N„C. 



B 



PT 



N C C A 



N r . 



C\ 



C 2 



B 



+ X eff M KK 



16n 2 K,J V 1Gtc 2 k,J \16tt 2 k, 
so that the eigenvalue of u> will be complex in a finite range of p when 



N 



(cl-c B ) + x 



l ff M 2 KK < 



(3.21) 



(3.22) 



167T 2 K / 

This determines the region of (Cb, Ca) having instability toward chiral magnetic spirals, 
which is a region defined by a hyperbola that is drawn in Figure HJ We will shortly confirm 
numerically that this rough picture indeed captures the true phase diagram qualitatively. 

It is important to notice that the purely baryonic chemical potential of Cb ^ with 
Ca = is stable against chiral magnetic spirals, and in fact we mathematically prove 
its stability in the Appendix. Therefore, previous Son-Stephanov chiral spiral phase is 
not affected by our results. This seems to indicate that axial chemical potential is the 
main cause of chiral magnetic spirals at strong coupling regime, while the weak coupling 
picture in Ref.[28] has both chemical potentials inducing chiral magnetic spirals. One 
possibility is that the dispersion relation of excitations that was used in Ref. J2B] might be 
drastically modified at strong coupling in the case of baryonic chemical potential, while 
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Pi P2 X p 

Unstable modes 



CO 




(a) Stable phase 



(b) Unstable phase 



(c) Phase boundary 



Figure 2: Spectral curves in (p, uj) in the case of (a) stable phase, (b) unstable phase, and 
(c) at the phase boundary. 

the case of axial chemical potential keeps its main features to have the phenomenon. We 
leave further understanding of this distinction to the future. 

To identify the true phase boundary of our interest in our numerical study (the hyper- 
bola in the case of the above approximation), the above rough analysis in fact provides 
us a very important clue how to search for the phase boundary numerically : consider 
the shape of dispersion relation of real uj and p as we cross the phase boundary. In the 
region without chiral magnetic spirals, that is, where the condition (I3.22p is not satisfied, 
the dispersion relation or spectral curve in real (p, u) space looks like a hyperbola whose 
transverse axis (meaning the line joining two focal points) is vertical. In real numerical 
situations, the true dispersion relation will be distorted and tilted, but the relevant feature 
is that the curve spans whole range of p while there is a gap in u. See Figure [2(a) for a 
schematic description. This is reasonable because we have real u for every p, and there is 
no instability for any p. In the opposite case of region with chiral magnetic spirals, that is, 
where (J3.22)) is satisfied, the dispersion relation curve would look like a hyperbola whose 
transverse axis is now horizontal, so that there is a gap in p instead while u is locally 
continuous. For p sitting between this gap, one has a complex eigenvalue for u, and hence 
instability to chiral magnetic spirals. See the Figure 12(b) . At the phase boundary, one 
expects something like the Figure E(c). What happens is that the two states with real 
eigenvalues of u in Figure E(a) become degenerate at the phase boundary, and then split 
into two states with complex eigenvalues of u conjugate to each other as one moves to 
the Figure 12(b). These states can no longer be seen in real {jp, oS) dispersion relation we 
are considering. Therefore, one can numerically distinguish the two phases by looking at 
dispersion relations in real (p,u) space as one changes parameters (Cb, Ca). 
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(a) Stable phase, Ca = 1.9 (b) Unstable phase, Ca = 2.1 (c) Phase boundary, Ca = 2.03 

Figure 3: Spectral curves of the lightest mode in the re-scaled (p, w) space with B = 1 
and Cb = 0. This confirms numerically the previous expectation in Figure El 



To illustrate how it works, we present an exemplar case of purely axial chemical poten- 
tial Ca 7^ with Cb = as in Figure EJ For convenience, we introduce the dimensionless 
variables 



P 



P 



UJ 



w 



M KK M KK 

so that the equation (13. 15[) reads as 



, [Cb,C/ 



%k 2 kM 2 kk 



{Cb, Ca) 



(3.23) 



r2\--. 



(W 



- f)A& + d z \{1 + Z 2 )d z A^] ± {pF 0Z + wF 3Z )A& = 



with 



Fqz 
F-sz 



1 



1 + Z 2 
1 



1 + Z 2 



Ca cosh(i? arctan Z) + Cb sinh(i? arctan Z) 



Cb cosh(S arctan Z) + Ca sinh(I? arctan Z) 



(3.24) 

(3.25) 
(3.26) 



Note that A* ' can be obtained from A^ by Ca — > —Ca and Cb -> —Cb- Also, there is 
a parity symmetry 



Z -»■ -Z 



Cn^C 



B 



Ca — > — C a 



P^ -P, 



w —7- w 



(3.27) 



so it is enough to study the equation of A^ in the first quadrant of {Cb, Ca) 1° identify 
the full phase diagram. Our numerical strategy to find the spectral curves in Figure 
E] is the 2D shooting method. We shoot from the initial point at Z^ ~ — oo with the 
normalizable boundary condition A^ + \Zi) = jr in ( 13.1 6p . and numerically integrate f )3.24p 
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Figure 4: Numerical phase diagram of chiral magnetic spiral at B = 1 and B = 3. It 
conforms to our expectation that a larger magnetic field should enhance the instability 
to chiral magnetic spiral. 



to Zf ~ +oo to find A( + '(Zf). Figure|3]are simple contour plots of A^ + \Zf) in the (p, w) 
space where lower values are darker. Only the right values of (p, w) will lead our shooting 



to be on vanishingly small values of 



Al+HZ 



fj 



at Zf ~ oo. From the Figure El^c) we can 
obtain one phase boundary point (Cb,Ca) = (0,2.03). By varying C B we search for the 
corresponding Ca showing the critical behavior such as Figure [3]^c), and this determines 
the phase boundary between chiral magnetic spiral phase and the stable homogeneous 
phase. 

The resulting phase diagrams of chiral magnetic spiral for two different values of 
magnetic field B = 1 and B = 3 are shown in Figure HI This is a numerical demonstration 
of the Figure HJ The red solid curve is for B = 1 and the blue dotted one is for B = 3. 
The larger magnetic field enhances the instability, so the region of chiral magnetic spirals 
gets bigger as expected. The full phase diagram is reflection-symmetric on Ca and Cb 
axis. This is a numerical corroboration of the chiral magnetic spiral in the strong coupling 
regime, at least for sufficiently large axial chemical potentials. 

To see more explicitly how complex eigenvalues w emerge, we trace them following the 
parameter space of Figure [3j i.e. B = 1, Cb = 0, from Ca = 1.9 to 2.1. As an example 
let us choose p — 0.85 which is a critical point shown in Figure [3]^c). Since all parameters 
except w are given, the eigenvalues w can be determined again by 2D shooting method, 
but now in (Re[w], Im[w;]) space. Some numerical results are plotted in Figure |5Ja). At 
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(a) Eigenvalues: numbers are Cas. 
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(b) Eigenfunction for Ca = 2.1, w = 0.19i. 



Figure 5: Eigenvalues and an explicit eigenfunction for B — 1, Cb = 0,p = 0.85. This is 
a concrete numerical proof of existence of instability toward chiral magnetic spiral. 

Ca = 1.9 there are two real eigenvalues at (w = ±0.26), which agree to the values read 
from Figure [3]^a). As Ca increases these two real eigenvalues approach to each other, 
and at Ca = 2.03 they merge at zero, as expected from Figure E](c). After passing this 
phase boundary point (Ca = 2.03), new pair of two imaginary eigenvalues appear and 
grow apart with increasing Ca- At Ca = 2.1 we find w = 0.1 9z of which complex- valued 
eigenfunction is shown in Figure |5^b). 

As a final comment, our analysis is based on the lowest mode spectral curve instability. 
There are towers of higher excited modes in the eigenvalue equation of ( 13.151) . and in 
principle there will be additional higher mode instability when they become unstable. 
Because these higher modes are more massive, the instability for them requires larger 
value of Ca than is needed for the lowest mode, and it happens always within the unstable 
region already determined by the lowest mode spectral curve. Thus these higher mode 
instability does not change the phase diagram. See Figure [6] for numerical illustrations 
on this point. 

4 Summary and future directions 

The fate of dense matter under extreme conditions is an interesting question relevant to 
high energy astrophysics and heavy-ion collisions. One such condition that we study is the 
application of large magnetic field, which indeed happens in the magnetars/neutron stars 
or off-center RHIC collisions. Especially, the presence of strong magnetic field provides a 
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Figure 6: Spectral curves in (p, w) at B = 1 and C# = 0, including a few higher modes too. 
After the lowest mode becomes unstable at Ca — 2.03, the next excited state becomes 
unstable at C A = 4.165. 

nice arena for physics of triangle anomaly to play important dynamical roles in the system. 
One of the most basic questions is the phase diagram of the favored ground state in terms 
of various conditions such as temperature, baryonic/axial chemical potentials, and the 
strength of magnetic field, etc. While some regions of this phase diagram are under 
control of perturbative QCD computations, there are other interesting regions whose 
characteristic scales are much lower than Aqcd, and strong interactions are crucial in 
understanding their properties. Holographic QCD based on the strong coupling expansion 
in large N c limit is one useful tool to study /gain insights on some of these physics, although 
there are several draw-backs and limitations too. 

We study the ground state of the system under a large magnetic field at zero temper- 
ature, confined/chiral-symmetry broken phase in the presence of both baryonic and axial 
chemical potentials. One possible phase that was proposed by Son-Stephanov [23j is a 
homogeneous distribution of axial/baryonic currents in the form of meson super-current 
along the direction of the magnetic field, whose magnitudes are determined exactly by 
triangle anomalies of QCD. On the other hand, Basar-Dunne-Kharzeev [28] recently pro- 
posed an alternative ground state, at least for sufficiently large magnetic field, where 
there are additional condensates of transverse currents that form spiral shapes along the 
magnetic field direction, and hence called chiral magnetic spirals. A basic assumption in 
their argument is the existence of fermionic quark quasi-particles in the effectively reduced 
model to 1+1 dimensions, which might be justified for sufficiently large magnetic field. 
However, for a more modest value of magnetic field, they are strongly coupled/confined, 
and it is not clear whether the picture remains effective in the strongly coupled regime. 
The purpose of this paper is to answer this question in the framework of holographic QCD 
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proposed by Sakai and Sugimoto. 

Our results indicate that the chiral magnetic spirals are indeed favored against the 
homogeneous phase at least for sufficiently large axial chemical potential. There is an issue 
here whether the relevant axial chemical potential that causes the phenomenon is within 
the validity regime of the model analysis. Especially, we used 2-derivative truncation of 
the DBI action on D8-branes neglecting higher derivative corrections. As our results in the 
previous section involve reasonable order 1 numbers in terms of dimensionless variables, 
we expect that our results do pass this issue safely. However, a more thorough analysis 
on this aspect, particularly interesting one being whether the same conclusion will be 
obtained if one uses the full DBI action instead, will be a valuable future direction. 

Our results do not have chiral magnetic spirals in the purely baryonic chemical po- 
tential, which in fact differs from the original proposal in Ref.[28] that the phenomenon 
should happen in both baryonic and axial cases. It would be interesting to understand 
why one case survives in the strong- coupling regime while the other does not. 
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A No chiral magnetic spirals for pure baryonic case 

In this appendix we will show rigorously that for pure baryonic case, we do not have any 
chiral magnetic spirals. Our starting point for this will be the master equation (" 13. 15[) 
replacing F oz and F 3Z with their purely baryonic expression ( (12. lip with Ca = 0). Let 
us first expand normalizable A ± ( y Z) in the complete basis of normalizable eigenf unctions 
f n {Z) satisfying the eigenvalue equation 

(1 + Z 2 )ld z \{l + Z 2 )d z f n (Z)} = -XlUZ) , (A.28) 



where < Xf < A 2 , ■ ■ ■ and the normalizable basis functions f n {Z) satisfy the completeness 
condition 

J(l + Z 2 )-^f: m (Z)fn(Z) = 5 nrn . (A.29) 

These are in fact wave- functions of vector mesons one can find in Refs.pTJ [33] . 
Because {f n (Z)} are complete, we expand 

A±{Z) = Y J <UZ) , (A.30) 

n>\ 

such that Xm>i l a n| 2 = I; so that A ± {Z) is normalized as 

f A±{Z)A ± {Z){1 + Z 2 )-^ = 1 . (A.31) 

Using the above equations and integrating the master equation in purely baryonic back- 
ground after sandwiching it with A ± (Z), we get the following 

{u 2 -p 2 )-M 2 KK b±{pa 1 -uja 2 ) = , (A.32) 

where b = J2 n >i (l a n| 2 A^) > Af > and 

N C C B f „-,,„* „+ skmfihan- 1 ^) 



- - i£/^W : 



1 + Z 2 
N C C B f ,-_,_,„, J+/r ^cosh( 

«2 



N Cn r - cosh fi? tan" 1 (Z)) 

^ A±(Z)AHZ) \ 79 ^ . (A.33) 

Stt^k j 1 + Z^ 



1 + Z 2 

Since cosh (Stan -1 (Z)) > sinh (i?tan _1 (Z)J, it is clear that \a 2 \ > \cti\- After complet- 
ing the squares in (1A.32[) . we get the following 

(« * f ) 2 = ( P T f )' + (^) + M^ (A.34) 

Since all the terms in the RHS of the above equation are positive, it is clear that u is 
always real in purely baryonic case and hence there is no instability and hence no signal 
of chiral magnetic spiral. 
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